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This article is a slightly modified version of the authors' contribution to Scattering, edited by Roy 
Pike and Pierre Sabatier, to be published by Academic Press. 



I. WHAT IS A BLACK HOLE? 

Black holes, objects so compact that not even light 
can escape their gravitational pull, are among the most 
intriguing concepts of modern science. As Kip Thornc 
wrote in 1974: "Of all the conceptions of the human 
mind from unicorns to gargoyles to the hydrogen bomb 
perhaps the most fantastic is the black hole: a hole in 
space with a definite edge over which anything can fall 
and nothing can escape; a hole that curves space and 
warps time." Evidence that these exotic objects exist 
in our universe has been mounting since the discovery 
of the archetypal black hole Cygnus XI in 1971. Today, 
the presence of supramassive black holes (several million 
times as massive as our sun) at the centre of many galax- 
ies, and smaller black holes (5-10 times as massive as the 
sun) in X-ray binary systems, is generally accepted. 

From the mathematical point of view, a black hole is 
a spacetime defined by a four-dimensional metric ten- 
sor, that is, a solution of Einstein's equations, with char- 
acteristic properties. The simplest black hole is spher- 
ically symmetric and non-rotating. It is known as the 
Schwarzschild black hole, and is described by the metric 
tensor g a b which is given by 



ds 2 = /J /J 9abdx a dx h = g a bdx a dx b = 

a=l 6=1 

A r 2 

= — -dt 2 + —dr 2 + r 2 d9 2 + r 2 sin 2 9 dip 2 (l.l) 

Here ds is the infinitesimal element of proper time, the 
coordinates are {x a } — (t,r,6,(p), and we have defined 
A = r 2 — 2Mr where M is a constant given by M = 
GMbh/c 2 where Mbh is the black hole mass and G is 
Newton's constant. In the following we will always choose 
so-called geometrised units in which both the speed of 
light c and G are equal to unity. We will also assume 
that repated indices indicate summation, as in 



yj) 



We are interested in scattering problems involving 
black holes. In such problems the curvature of spacetime 
enters not only at the level of the boundary conditions, 
but also in the equations describing the propagation of 



the various wave-fields (scalar, electromagnetic or grav- 
itational) that we may be interested in. Therefore the 
problem of scattering by black holes has more in com- 
mon with scattering in media with a non-constant index 
of refraction than scattering by a physical object. In the 
case of a black hole, it is the curvature of space-time itself 
which is doing the scattering. 

In this article we aim to summarise work in this re- 
search area, and relate them to results in familiar con- 
texts such as quantum scattering. We will discuss pos- 
sible diffraction effects, resonances, and digress on some 
of the peculiarities of the black-hole problem. It is useful 
to begin our discussion with the simplest of scattering 
problems involving black holes; namely, the bending of 
light rays propagating in a black hole spacetime. 



II. CLASSICAL TRAJECTORIES 



Photon trajectories outside Schwarzschild black 
holes 



Let us consider an astronaut who, while piloting his 
spaceship towards a Schwarzschild black hole (Figure |l|), 
shines a laser directly out his window, in the positive tp 
direction. The trajectory of the laser beam can be found 
by solving the equation for a null geodesic (a line that is 
"as straight as possible" in the curved spacetime) . Hence 
we are looking for a solution to the geodesic equation 
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(2.1) 



where the Christoffel symbols are functions of the 



coordinates {x a }. In order to be null, our geodesic must 
satisfy x^x^ = 0. Now , if w e use the Schwarzschild coor- 



dinates introduced in (1.1) we can cast the equation for 
a null geodesic into the following form: 



dip = _t_ 

dr r 2 
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(2.2) 



where b is the impact parameter defined by b = L/E 
where L and E are the angular momentum and the en- 



2 



ergy associated with the photon. From equation ( |2.2| ) we 
can deduce the properties of various light trajectories. 



FIG. 1: An astronaut maneuvers his rocket ship near a non- 
rotating black hole. 




As the astronaut nears the black hole, the laser beam 
is deflected more and more by the spacetime curvature — 
see Figure 0. The curves in Figure ||| are the solutions of 
equation (2.2) with cp' (b) = and b ranging from 2.5M 
to 5M. It should be noted that r = 3M corresponds 
to a circle. This is known as the unstable photon orbit, 
and its existence means that our astronauts laser beam 
will circle the black hole and illuminate his neck! After 
passing r = 3M the astronaut finds that the laser beam 
is always deflected into the black hole. 



light at various angles. He then finds that there is still 
a (small) range of angles at which the light beam can 
escape the black hole (see Figure 0). 



FIG. 3: More light trajectories in the Schwarzschild space- 
time. At r = 2.1M, the astronaut finds that a light beam 
from his spaceship can escape the black hole for only a small 
range of angles. As r — > 2M, this becomes a single point. 




Once the astronaut has resumed his fall towards the 
black hole, and reached r — ► 2M, the solid angle into 
which he must shine his laser in order for the light to 
escape the black hole has shrunk to a single point. He 
must aim it in the positive r direction. However, even 
if he hovers at constant r = 2M, the ratio of the light 
frequency received (y') by his home planet at (say) r' — 
oo to the frequency emitted (v) is the ratio of the proper 
times at the two points. This follows from the formula 
for the gravitational redshift: 



FIG. 2: Light trajectories in the Schwarzschild geometry for 
various values of the impact parameter b. This shows what 
happens to a laser beam which shines out the window of a 
spaceship as it plunges into a (Schwarzschild) black hole. Note 
the circular 'photon orbit' at r = 3M. The (grey) circle at 
r — 2M represents the event horizon of the black hole. 



v' = 1/(1 - 2M/r). (2.3) 

From this we see that when r — > 2M the light is com- 
pletely redshifted away. Therefore any light emitted after 
the astronaut reaches this limit, whatever the direction, 
remains inside the black hole. Our space traveller has 
reached the so-called event horizon, and he can neither 
escape the black hole nor alert a rescue team of his fate. 




r/M 



When he reaches r = 2.1M, the astronaut applies his 
rockets in such a way that the spaceship hovers at a con- 
stant distance from the black hole, and tries to shine his 



B. Bending of Starlight 

The above discussion illustrates some of the extreme 
effects that the curvature of a black-hole spacetime might 
have on light trajectories. Still, the ideas are relevant also 
in a more familiar setting. In fact, the first experimental 
verification of Einstein's theory of general relativity was 
the measurement of the bending of starlight by the grav- 
itational field of the sun during a solar eclipse of 1919. 
While the sun is certainly not a black hole, the metric 
tensor exterior to the surface of the sun is accurately de- 
scribed by ([□]). Thus, in a sense, the first test of general 
relativity was also the first 'black-hole scattering' exper- 
iment. 
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FIG. 4: The earliest 'black-hole scattering' experiment: de- 
flection of starlight by the sun. 
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Let us assume that M/r is small (<< 1) but not neg- 
ligible. If we introduce a new variable 



1 



y 



i 



M 
r 



we can approximate our equation ( p.2[ ) as 
1 + 2My 



dy 



0{M 2 /r 2 ^ 



(2.4) 



(2.5) 



Integration of this leads to 



2M -t / , ( 1 , 

ip - <p = — + sin \by) -2Mi—-y 2 



1/2 



(2.6) 



The initial trajectory was such that y — > 0, which means 
that tp — > (po, defining the incoming direction. The small- 
est value of y that the light ray will ever reach corre- 
sponds to y = 1/b when ip — ipo = 2M/b + ir/2. By 
symmetry the ray will be deflected by an equal amount 
as it recedes to infinity. Hence, the final result will be 
tp — ifQ = AM/b + n, and since the result would have been 
ip — <y5o = ft if the light ray had followed a straight path, 
we have a total deflection of 



AM 



(2.7) 



This result was first calculated by Einstein, and as noted 
above, this deflection was measured by Eddington and 
collaborators in 1919. 



C. Gravitational lensing 

With enhanced observational capabilities, ranging 
from large radio telescopes to the repaired Hubble space 
telescope, astronomers have found plenty of evidence that 
light often bends as it travels through space. This is 



known as gravitational lensing. Before discussing some 
of these observations it is instructive to describe the sim- 
plest lensing situation; the case when light from a distant 
quasar is lensed by a localised mass distribution (treated 
as a point source with mass M) on its way to the ob- 
servatory. The relevant geometry is shown in Figure [|. 
We denote the true angular separation from the lens to 
be (3, consider a light ray that passes the lens at a min- 
imum distance £ (essentially the impact parameter from 
the previous section) and is deflected an angle a — 4M/£ 
by the spacetime geometry. Assuming that £ > > 2M we 
do not need to worry about diffraction effects, and can 
analyse the problem using simple geometry. Disregard- 
ing complicating factors one can show that the condition 
that the light ray reaches the observer leads to 



(3 = 6- 

where we have defined 

0- 



(Is 



D, 



(2. 



4M- 



Da 



D d D s 



and the distances -Dds, and D s are as shown in Fig- 
ure |^. 

FIG. 5: A schematic description of the simplest gravitational 
lensing geometry. Light rays are bent as the pass by a point 
source located at point M on their way from the distant source 
at S to the observer O. 




In other words, we solve 

9 2 - (39 ~ al = 

and get 



(3± 



Aa 2 



(2.9) 



(2.10) 



That is, we always get two solutions of opposite sign. 
This means that there will typically be one image on each 
side of the lens. A special case worthy of notice arises 
when the source, lens and observer are all aligned. Then 
we have (3 = and it is easy to realise that there is then 
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no preferred plane for the light rays to travel in. Thus the 
whole ring of angular radius \8\ — ao is a solution to the 
simple lensing equation. This is commonly known as an 
'Einstein ring' and, as is easy to see, such images can only 
occur in lensing by axially symmetric mass distributions. 

The above description is obviously idealised in many 
ways, and in analysing observed gravitational lenses one 
must consider much more complicated mass distribu- 
tions, as well as use detailed cosmological models. Still, 
the above example illustrates the basic principles, and 
we now turn to some actual observations of gravitational 
lensing. 

The first lensing candidate was observed using a radio 
telescope at Jodrell Bank in 1979 and is catalogued as 
0957+561. It is a typical example of a double image of a 
distant quasar. When the spectra of the two images were 
studied it was found that they were remarkably similar, 
but redshifted to slightly different extent. It was con- 
cluded that the two images were extremely unlikely to 
correspond to separate individual quasars. 

Since this first discovery, many other double image 
systems have been observed. More complicated sys- 
tems, comprising of further images, have also been found. 
Among them are the triple image lenses 2016+112 and 
0023+171. Perhaps the most unique multiple image case 
found so far is 2237+0305, more commonly known as 
the "Einstein cross" . In this beautiful system, shown in 
figure ^, one can see four distinct images. 



FIG. 6: The famous "Einstein Cross" is a case where a grav- 
itational lensing by a massive galaxy (central image) leads to 
four distinct images of a distant quasar. This picture was 
taken by the William Hershel Telescope in August 1994. 




complete rings. One interesting example of the latter is 
MG1131+0456. 

As our observational capabilities continue to improve 
the list of lensed systems is rapidly growing. New obser- 
vations provide challenges for the theorists that want to 
deduce the geometry of the lensing mass distribution as 
well as understand the nature of the original light source. 
Ideally, one would also like to be able to use lensing ob- 
servations of distant quasars to also deduce information 
about cosmology. Considerable progress in these direc- 
tions have been made in recent years, and our under- 
standing should continue to improve with the observa- 
tional data. For further details, we refer the reader to 
the monograph by Schneider, Ehlers and Falco Q. 



III. WAVE SCATTERING 

Having discussed some classic examples and exciting 
observations of the scattering of light by massive bodies 
we now turn to the issue of possible diffraction effects. 
Essentially, most models of the gravitational lensing phe- 
nomenon are based on "geometrical optics" . Therefore 
one would not expect these calculations to yield any in- 
sight into possible wave phenomena. However, as we will 
see, the extreme nature of black holes lead to the exis- 
tence of many complicated diffraction effects. To under- 
stand these it is essential that we develop a framework for 
studying the scattering of waves by black holes. From an 
observational point of view, our main interest will be fo- 
cussed on electromagnetic and gravitational waves. The 
latter are particularly interesting since a new generation 
of gravitational-wave detectors is due to come online in 
the next few years. It is generally believed that these will 
make the long-heralded field of gravitational- wave astron- 
omy a reality, and that they will allow us ot make detailed 
observations of the physics in the immediate vicinity of 
a black hole ||. 

In order to introduce the various concepts involved in 
studies of the scattering of waves by black holes we will 
consider the relatively simple case of scalar waves. This 
may seem like a peculiar choice given that no massless 
scalar fields have yet been observed in nature. However, 
it turns out that the main equations governing a weak 
electromagnetic field, or gravitational waves, in a curved 
spacetime are essentially the same as the sc alar field wave 
equation (see for example equation (8.18)). Hence, the 
scalar field serves as a useful model. See Chandrasekhar 
S for an exhaustive study. 



The above examples are all likely examples of lens- 
ing by 'point sources'. For extended mass distributions, 
like galaxies, one would expect to see also arcs and in 
some unique cases almost complete Einstein rings. The 
first arcs (Abell 370 and C1224-02) were actually found 
when gravitational lenses were still considered a mere 
theoretical possibility. Since then many further exam- 
ples of lensed arcs have been discovered, as well as nearly 



A. Scalar Fields in the Schwarzschild geometry 

Let us consider a scalar field $ propagating in the 
Schwarzschild spacetime, as described by (1.1). The 
equation governing the evolution of the scalar field is 
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where g is the determinant of g^ v . Since the 
Schwarzschild spacetimc is spherically symmetric, we 
may assume a Fourier decomposition 



1 OO I 



<!<7"(«, y )0,„,(u,r) (3-2) 



i=0 ' 



where YJ m are the spherical harmonics. If the boundary 
conditions are cylindrically symmetric, as they would be 
in the plane-wave scattering problem, the ip — axis (pr- 
axis) can always be chosen to be the axis of symmetry. 
We can therefore assume that <pi m to be independent of 
(p and perform the sum over m, writing 



(3.3) 



where P; is a Legendre function. 

In the Schwarzschild spacetime the wave equation for 
the scalar field reduces to the following Schrodinger-type 
equation for </>;: 



[u; 2 - V{r)} = 0, 



(3.4) 



where the so-called tortoise coordinate is defined by 

r - 2M d 



d 
dr* 



This integrates to 



dr 



2Mlog(— - 1 
\ 2A1 



(3.5) 



(3.6) 



and we see that introducing the tortoise coordinate cor- 
responds to "pushing the event horizon of the black hole 
away to — oo". The effective potential is explicitly given 
by 



V(r) 



2M 



1(1 + 1) 



2M 



(3.7) 



It is positive definite and has a single peak in the range 
r* e [— oo, oo], see Figure 0. 

A black hole is distinguished by the fact that no infor- 
mation can escape from withi n the event horizon. Hence, 
any physical solution to (3.4) must be purely ingoing at 
the event horizon, that is, at r = 2M (r* — > — oo). There- 
fore we seek solutions to (3A) of form (for a given fre- 
quency oj) 



oo 
— oo 
(3.8) 

where the amplitudes of the scattered and the transmit- 
ted waves, Si and TJ, remain to be determined. Clearly, 
problems involving waves scattered from a Schwarzschild 
black hole share many features with scattering problems 



FIG. 7: A schematic description of the scattering of waves 
in the Schwarzschild background. The effective potential of 
equation (3.4) is shown as a function of r*. The event horizon 
of the black hole is located at r* = — oo. An incident wave I is 
decomposed into a transmitted component T and a scattered 
component S. 
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in quantum theory. Hence, we can adopt standard tech- 
niques to evaluate Si and T/. By conservation of flux it 
follows that 



\Ti\' = l-\S, 



(3.9) 



Hence, we need only determine either Si or T/. Typical 
results for Si are shown in Figure |^. 

The nature of Si can be understood from the follow- 
ing observations. For u> « 2M the wavelength of the 
infalling wave is so large that the wave is essentially un- 
affected by the presence of the black hole. It is only if 
we "aim" the wave straight at the black hole (recall that 
the impact parameter follows from b = L/E ~ l/u) that 
we can get an appreciable effect. Hence, we expect to 
have Si — > 1 as uj — > 0. For large frequencies oj » 2M, 
the situation is the opposite and we expect to find that 
Si — > as uj — ► oo. Thus, high frequency waves will be 
absorbed unless they are aimed away from the black hole. 



FIG. 8: Scalar wave scattering coefficient Si for I = 0,1, 
as a function of 2Muo. 
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In studies of dynamical black holes one often needs 



to construct the general solution to (3.4). One can do 
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this using two linearly independent solutions. These are 
customarily normalised in a way that differs slightly from 
(3.8). A first solution (essentially (|3.8[ )) is such that the 
amplitude of the waves falling across the event horizon 
are normalised to unity, and one requires the amplitudes 
of out- and ingoing waves at infinity. This solution is 
sometimes called the "in" -mode and it can be written 



A out (u)e iur - + A- m (u)e 



— oo 
+00 



(3.10) 

Given this solution, a second linearly independent one 
corresponds to waves of unit amplitude reaching spatial 
infinity. This is the "up" mode, and it follows from 



up f B out (co)e^ r ' + B- m {u)e 



1 

r* 



-00 
-00 



(3.11) 

The nature of these two solutions is illustrated in Fig- 
ure H 

FIG. 9: The nature of two linearly independent solutions to 
the scalar wave equation outside a black hole. The in-mode 
corresponds to purely ingoing waves crossing the event hori- 
zon {H + ), while the up-mode corresponds to purely outgoing 
waves at spatial infinity 




the black hole. The various arrows represent the path 
followed by the incident, transmitted and reflected wave 
fronts. 



B. Plane wave scattering 

In order to better understand the physics of black holes 
we want to formulate a scattering problem analogous to 
that used to probe the nature of (say) nuclear particles. 
We want to let a plane wave fall onto the black hole, 
and investigate how the black hole manifests itself in the 
scattered wave. This is obviously a model problem, given 
that we cannot expect to ever be able to compare the 
calculated scattering cross sections to real observations, 
but it is still instructive. In particular, it will lead to an 
unveiling of a deep analogy between black hole physics 
and well known phenomena such as glory scattering. 

However, in formulating this problem we immediately 
face difficulties. What exactly do we mean by a plane 
wave in a curved spacetime? It turns out that we can 
answer this question by appealing to the analogous prob- 
lem of Coulomb scattering. As with the charge in the 
Coulomb problem, the black hole contributes a long- 
range potential that falls off as 1/r at large distances. 
The effect of such a long-range potential on the "plane" 
wave can be accounted for by a simple modification of 
the standard (flat space) expressions for the scattering 
amplitude. In the black-hole case we essentially need to 
introduce r* in the phase of the plane wave, and we ob- 
tain 
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-00 in the case of scalar waves. 



In 



(3.12) 



For those unfamiliar with conformal diagrams, a short 
explanation of Figure is necessary. Suppressing the an- 
gular coordinates [9, <p), we make a coordinate transfor- 
mation in (t, r) such that the half-infinite space exterior 
to the black hole is mapped onto a finite portion of the 
plane. The transformation is chosen so that the light 
cones always intersect at an angle of 45° — much like 
Mercator's projection of the earth distorts the shape of 
the continents but preserves the directions North-South 
and East- West. The four 'points at infinity' are mapped 
into the diagonal edges as follows: The diamond-shaped 
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figures represent the whole of the space-time exterior to 



C. Phase-shifts and deflection function 

Having defined a suitable plane wa ve, t he scattering 
problem involves finding a solution to (3.4) , i.e., identi- 
fying the asymptotic amplitudes A- ln and A out for a given 
frequency lu. Then we can extract the scattered wave by 
discarding the part of the solution that corresponds to 
the original plane wave. The physical information we are 
interested in is contained within the scattering amplitude 
/(#), which follows from 



$ ~ $ P ianc + J -^-e iu}r ' , asr-,^ +00 



(3.13) 



Now letting 



$ _ $ 



plane 



2icur 



5^(2Z + 1) [e 2U ' -l]P;(cos0) 
+00 , (3.14) 
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define the (complex-valued) phase-shifts 61 it is straight- 
forward to show that 



Si = (-i) 



1+1 A)ut 



(3.15) 



From this it follows that the scattering amplitude, that 
contains all the physical information, is given by 

oo 

= Y~ ^ {21 + 1} ^ - 1] P ' (cos e) ■ (3 ' 16) 

1=0 

When discussing the physical quantities that follow 
from a set of phase-shifts it is natural to use Ford and 
Wheeler's excellent description of semiclassical scatter- 
ing from 1959 [Q. In the semiclassical picture the impact 
parameter 6 is given by 



b=[l 



1\ 1 



21 u> 



(3.17) 



Then each partial wave is considered as impinging on the 
black hole from an initial distance b away from the axis. 

In this description much physical information can be 
extracted from the so-called deflection function 6(J). It 
corresponds to the angle by which a certain partial wave 
is scattered by the black hole, and is related to the real 
part of the phase-shifts; 



9(0 = 2-Re Si 
al 



(3.18) 



Here / is allowed to assume continuous real values. We 
can obtain approximations for 6(Z) in some limiting 
cases. For large values of the impact parameter, 6, one 
would expect the value of the deflection function to agree 
with Einstein's classic result 9 « — 4M/6, that we de- 
rived earlier. As can be seen in Figure [to] this is cer- 
tainly the case. A numerically determined rapidly 
approaches the approximate result as / increases. 

A second approximation is intimately related to the ex- 
istence of a glory in black-hole scattering. Whenever the 
classical cross section diverges in either the forward or 
the backward direction a diffraction phenomenon called 
a glory arises. This phenomenon is well-known in both 
optics and quantum scattering. In general, backward glo- 
ries can occur if 6 < — tt for some values of b. Whenever 
the deflection function passes through zero or a multiple 
of 7r we have a glory. In the black-hole case one would ex- 
pect glory scattering to be associated with the unstable 
photon orbit at r = 3M. This essentially means that we 
would expect a logarithmic singularity in the deflection 
function to be associated with the critical impact param- 
eter b c = 3\/3M . This feature is obvious in Figure [Tc| . 
Many years ago Darwin || deduced an approximate re- 
lation between the impact parameter and the deflection 
function close to this singularity; 

6(6) w 3\/3M + 3.48Me" e . 



(3.19) 

If we invert this formula and use ( 3.17|) we get 6 as a 
function of /. As can be seen in Figure [10| this approxima- 
tion is in excellent agreement with the deflection function 
obtained from the approximate phase-shifts. 



FIG. 10: The deflection function O (solid) is shown as a 
function of I for u>M = 10. For large impact parameters (large 
I) the approximate results approach the Einstein deflection 
angle — 4M/6 (dashed). A logarithmic singularity in O is 
apparent at the critical impact parameter (l c ~ 51.5 here). 
This feature is associated with the unstable photon orbit at 
r = 3M. Also shown (as a dashed curve) is an approximation 
obtained by inverting Darwin's formula. 
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D. The black- hole glory 

We now want to proceed to calculate the s cattering 
amplitude through the partial-wave sum (3.16). In do- 
ing this, we must proceed with caution since, as in the 
Coulomb problem, the sum is divergent. This problem is 
normally avoided by introducing a cutoff where the re- 
mainder of the true partial- wave sum is replaced by ana- 
lytic results for a limiting case. In essence, we extr act t he 
contribution from large impact parameters from (3.16) , 
i.e., replace it by 



f(9) = f N (6) + f D (9) , 



(3.20) 



where the long-range (Newtonian) contribution is given 
by 



MB) 



r(i 

M- v 



2iMco) 



T(l + 2iMuj) 



sin — 
2 



-2+4iMui 



(3.21) 



and 



f D (0) 



2iuo ' 

2=0 



^(2/+l) 



2lS, 



,2Mf 



P; (cos 6»), (3.22) 



is the part of the scattering amplitude that gives rise to 
diffraction effects. The Newtonian phase-shifts 5^ follow 
from (cf. the standard Coulomb expression) 



o -IN 



I\Z + 1 - 2iMu) 



(3.23) 



r(Z + 1 + 2iMuj) 

The sum in ( [3.22 ) is convergent, and we can readily de- 
termine the desired physical quantities. The quantity of 
main physical interest is the differential cross section - 
the "intensity" of waves scattered into a certain solid an- 
gle. It follows from the well-known relation 



da 



um 2 



(3.24) 
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A typical example of a scalar wave cross section is shown 
in Figure [□]. This provides a beautiful example of the 
black-hole glory (the regular oscillations at large angles) . 
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FIG. 11: The scalar wave differential cross section for ujM = 
2.0. The oscillations on the backward direction corresponds 
to the black-hole glory. Also shown is an ap proxi mation to 
the glory oscillations following from equation (3.2E). 



It is appropriate to discuss the details of the cross sec- 
tion shown in Figure [ll] in somewhat more detail. First 
of all, one would expect the long-range attraction of the 
gravitational interaction to give rise to a divergent focus- 
ing (~ 6~ ) in the forward direction. Secondly, we expect 
to find that interference between partial waves associated 
with the unstable photon orbit, i.e., with impact param- 
eters b w b c — 3y/3M, gives rise to a glory effect in the 
backward direction (cf. Figure |l^). This effect is, of 



course, prominent in Figure 11. The glory oscillations in 
the cross section can be approximated by0 



da 
dtt 



= 2TTU>b 2 



Jq (uib sin ( 



(3.25) 



When combined with the Darwin formula (3.19) this pro- 
vides a good approximation whenever toM » 1 and 
\9 - 7r| << 1, see Figure [ll]. 

According to the predictions of geometrical optics one 
might expect to find glory oscillations not only in the 
backward direction, but also in the forward direction. 
Partial waves associated with the critical impact param- 
eter may be deflected any multiple of 7r and so give rise to 
diffraction close to both (9 = and ir. Moreover, in Fig- 
ure [l^ we see that the deflection function passes through 

= for a value of I lower than that associated with the 
unstable photon orbit (l c ). This means that there should 
be also a forward glory. However, this feature drowns in 
the divergence of the cross section that is due to the large 

1 partial waves. Moreover, partial waves corresponding 
to I < l c are to a large extent absorbed by the black hole. 
Hence, the forward glory is exceptionally faint. 



IV. TIME-DEPENDENT SCATTERING — 
RESONANCES 

As we have already pointed out, the scattering theory 
of monochromatic plane waves by a black hole, though 
interesting, is not likely to have any experimental con- 
firmation soon. The only possibility (as we can see it) 
would be via diffraction effects observed in gravitational 
lenses. In order to discuss potentially observable scat- 
tering effects we must turn to the closely related prob- 
lem of the response of the black hole to an initial per- 
turbation. The perturbation response, the black hole's 
"fingerprints", will only be observable via gravitational 
waves. Such observations are still outstanding, but the 
should become reality in the next few years when large 
scale interferometric gravitational wave detectors (LIGO, 
VIRGO, GEO600, TAMA300) come online. A typical as- 
trophysical scenario that would produce copious amounts 
of gravitational waves is a black hole in the act of swal- 
lowing a neutron star at the endpoint of binary evo- 
lution The perturbation of the black- hole gravita- 
tional field would produce outgoing gravitational radia- 
tion which can be detected far from the black hole. We 
shall see below that this radiation has a characteristic 
spectrum which indicates the presence of the hole and 
can be used to deduce its mass and angular momentum. 
In other words, gravitational- wave observations will bring 
direct observations of the most elusive of our universe's 
inhabitants, the black holes. 



A. A Green's function approach 

We are interested in modelling the response of a black 
hole to a prescribed initial perturbation. As a model 
problem we continue with the scalar field, and consider 
the evolution of the field as an initial-value problem. 
That is, instead of prescribing the asymptotic charac- 
ter of the field (as in the plane-wave case) we suppose 
that we are given a specified scalar field at time t = 
and that we want to calculate the future evolution of the 
field. The time-evolution of a wave-field $i(r»,i) follows 
from 



+ / d t G(n,y,t)^ l (y,0)dy , (4.1) 



for t > 0. The (retarded) Green's function is defined by 
d 2 d 2 



dr 2 dt 2 



V(r) 



G{n,y,t) = 5{t)8{n-y) , (4.2) 



together with the condition G(r*,y, t) = for t < and 
the appropriate space boundary conditions. These con- 
ditions follow from 



dG 



+ iujG = , r 



2M 



(4.3) 
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dG 

— iuiG = , r — > oo . 



(4.4) 



Let us take G to be the one-sided Fourier transform of 



G: 



G(r*,y,t)e luJt dt 



(4.5) 



This transform is well defined as long as Im u> > 0, and 
the corresponding inversion formula is 



1 

2^ 



+ 00 + 2C 



G{r*,y,uj)e- lut du , (4.6) 



— OO + 'iC 



where c is some positive number (see Figure |l2|). The 
Green's function G(r*,y, cj) can now be expressed in 
terms of two linearly independent solutions to the homo- 
geneous equ ation (3.4). The two required solutions are 
( 3.10 ) and ( 3.11 ) as defined earlier, giving the Green's 
function: 



G(r*,y,u) = - — 



1 



2iujAi n (u>) 



<f>\ n (y,u)<j>1 P (r*,uj) , r„ > y 
(4.7) 



Here we have used the Wronskian relation 



W(lu) = - 4>7^- = 2iuM in (w) ■ (4. 



B. Quasinormal modes 

The initial- value problem can now, in prin ciple , be ap- 
proached by direct numerical integration of (3.4) for (al- 
most) real values of u> and subsequent inversion of (4.6). 
It has proved useful, however, to deform the contour of 
integration in the complex ui plane using Cauchy's theo- 
rem and rewrite the integral as a sum over residues plus 
a remainder integral. Our first task in this process is to 
find the position of the poles of the Green's func tion. 

The poles in the radial Green's function (4/7) are lo- 
cated at the simple zeros of Ai ni which we will denote 
{ujq\. When ui — u) q , the distinction between <p™ and 

disappears, so that, if we demand that our solution 
to (3.4) be both ingoing at the horizon and outgoing at 
infinity, and then solve the resulting relation for ui, we 
know that we have found a pole. 

It turns out that the corresponding frequencies, known 
as the quasinormal modes of the black hole (quasinormal 
because they are damped as radiation dissipates to infin- 
ity and across the event horizon), play a dominant role 
in the evolution of black hole perturbations. The first 
indication of this was found by Vishveshwara |8| . He re- 
alized that one might be able to observe a solitary black 
hole by scattering of radiation, provided the black hole 
left its fingerprint on the scattered wave. So he started 
"pelting" the black hole with Gaussian wave packets. By 



FIG. 12: The integral contour of equation ( |4.6| ). If t > r„, the 
contour can be deformed and the integral can be rewritten as 
a sum over the poles (crosses) and a remainder integral over 
the branch cut (zig-zag line). 



tuning the width of the impinging Gaussian Vishvesh- 
wara found that the black hole responded by ringing in 
a very characteristic decaying mode; the slowest damped 
of the black holes quasinormal modes. Subsequent work 
(in particular in numerical relativity) has shown that the 
quasinormal modes always play a prominent role in the 
dynamical response of a black hole to external pertur- 
bations (see Figure |l3|). Impressive results for head-on 
collisions of two black holes lead to signals that are al- 
most entirely due to quasinormal mode ringing. 

The actual determination of quasinormal mode fre- 
quencies is a far from trivial cal culation. The quasinor- 
mal modes are solutions to (3.4) that satisfy the causal 
condition of purely ingoing waves crossing the event hori- 
zon, while at the same time behaving as purely outgo- 
ing waves reaching spatial infinity. Assuming a tim e- 
dependen ce e" 1 "*, a general causal solution to (3.4) is 
given by ( 3.1C ) and a quasinormal mode corresponds to 
A- m = 0. To identify a mode-solution we must therefore 



be able to determine a solution that behaves as e lujr * as 
r* — > oo, with no admixture of ingoing waves. Assum- 
ing that the black hole is stable (we can in fact prove 
that this must be the case), no unstable mode-solutions 
should exist so we must require that a mode is damped 
according to an observer at a fixed locati on. T his means 
that Im uig < 0. The general solution (3.1C) is then a 
mixture of exponentially growing and dying terms. We 
must, out of all solutions, identify the unique one for 
which the coefficient of the exponentially dying solution 
is zero. Several methods have been devised to deal with 
this difficulty accurately ||. These methods have been 
used to investigate the entire spectrum for non-rotating 
black holes, and also to map out the behaviour of the 
first ten modes or so as the black hole spins up. The 
spectrum of gravitational- wave modes of a Schwarzschild 
black hole is shown in Figure [l4|. 

It is worthwhile to outline one of the most reliable 
methods (due to Leaver jl(|) for calculating black hole 
quasinormal modes. Write the desired solution to equa- 
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FIG. 13: A recreation of Vishveshwara's classic scattering 
experiment: The response of a Schwarzschild black hole as a 
Gaussian wavepacket of scalar waves impinges upon it. The 
first bump (at t = 50Af ) is the initial Gaussian passing by 
the observer on its way towards the black hole. Quasinormal- 
mode ringing clearly dominates the signal after t » 150M. 
At very late times (after t ~ 300A/) the signal is dominated 
by a power-law fall-off with time. This late time tail arises 
because of backscattering off of the weak potential in the far 
zone. As such, it is not an effect exclusive to black holes. A 
similiar tail will be present also for perturbed stars. 



tion (3.4) as an infinite sum 



4>T = (r-2M) p (2M/ry p e 



2 P „-p(r-2M)/2M 



n=0 



r-2M\' 



(4.9) 

where p = —i2Mu>. The recurrence relations between the 
a n are given by Leaver: 



a n a n+1 + (3 n a n + 7„a„_i = 0, 



(4.10) 



where 



a n = n 2 + 2n(p + 1) + 2p + 1 

f3 n = -[2n 2 + 2n(4p+l) + 8p 2 +Ap + 1(1 + 1) - s 2 + 1] 
Jn = n 2 + 4np + 4p 2 - s 2 

and s is the spin of the field. Now one can note that the 
coefficient A in has a zero whenever the sum ^ a n con- 
verges. This requirement translates into an continued- 
fraction equation involving the coefficients a, ft and 7: 



Here q — 1,2,3,... is the mode number. Solving this 
equation numerically (still a non-trivial task!) for u q 
gives the quasinormal modes (see Figure |l4| ). 



FIG. 14: The complex quasinormal mode frequencies corre- 
sponding to gravitational perturbations (for I — 2 and 3) of a 
Schwarzschild black hole. These correspond to the p ositi ons 
of the poles of the radial Green's function (equation 4.7) in 
the complex 2Mu plane. 
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C. Mode excitation 

Having located the quasinormal modes we want to 
evaluate the contribution of each mode to the emerging 
signal. Ideally one would like to be able to quantitatively 
account for the contribution to a signal from each indi- 
vidual quasinormal mode. Thus we want to con struct the 
mode-contribution to the Green's function (4.6), combine 
it with the releva nt in itial data an extract the correspond- 
ing signal using ( |4.l| ). To do this is largely a (rather in- 
volved) numerical exercise. In the end one finds that the 
quasinormal modes account for the main part of the sig- 
nal after a certain time, essentially the time it takes for a 
part of the initial data to travel from its original position 
y to the black hole, and then for the scattered wave to 
reach the observer at r*. Thus we expect the modes to 
generally dominate for (roughly) t — r* — y > 0. 

It is helpful to introduce a simplifying approximation 
at this stage. Let us assume that the initial data has 
support only far away from the black hole, and that the 
observer is also located in the far zone. Then we can 
can replace the solutions 



bi in (4/7) by their asymptotic 
behaviour at large r* and readily evaluate the mode con- 
tribution to the Green's function. Since A- ln (ui) has a 
simple zero at u> = u> q , it is useful to define a quantity a q 
by 



A in (u>) w (w - uj q )a q , 



(4.12) 



in the vicinity of the pole. Then it follows from the 
residue theorem (and the fact that modes in the third 
and fourth quadrant are in one-to-one correspondence, 
see Figure 113) that the total contribution from the modes 
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to the time-domain Green's function can be written 



G Q (r*,y,t) =Rc 



Here we have defined 



B„ 



q=0 



-iuj q (t-r,-y) 



(4.13) 



(4.14) 



The sum in (4.13) is over all quasinormal modes in the 
fourth quadrant of the w-plane. That this expression pro- 
vides an accurate representation of the mode-excitation 
(as long as our assumptions arc valid) has been demon- 
strated. 

At this point it is relevant to comment on the fact that 
the quasinormal modes are, even though there is an in- 
finite set of modes for each I, not complete. That is, a 
mode sum such as ( 4.13 ) should not be expected to repre- 
sent the entire black-hole signal for given initial data. It 
will typically not be useful at early times, and it cannot 
represent the power-law tail that dominates at very late 
times, see Figure [13L However, the mode-contribution 



is still highly relevant, and results like ( 4.15 ) can help 
us understand the dy namic s of black holes better. Fur- 
thermore, expression ( 4.13 ) allows us to study the con- 
vergence of the mode-sum in a simple way. It has been 
shown that (again under the assumptions of the "asymp- 
totic approximation" ) that the mode sum converges for 
t — r* — y > 0. 



which is accurate to within 10 percent for the fundamen- 
tal mode. 

Let us translate the results for the fundamental 
gravitational-wave quasinormal mode of a nonrotating 
black hole into more familiar units. We then get a fre- 
quency 



/ „ 12kHz ( ^ 
1 M 



(4.17) 



where Mq represents the mass of the Sun, while the as- 
sociated e-folding time is 



0.05 ms 



M 



(4.18) 



The quasinormal modes of a black hole are clearly very 
shortlived. In fact, we can compare a black hole to other 
resonant systems in nature by considering the quality 
factor 



Q 



Re oj a 



Im io 



(4.19) 



Our quasinormal-mode approximations then lead to Q ~ 
/. This should be compared to the result for the fluid 
pulsations of a neutron star: Q ~ 1000, or the typical 
value for an atom: Q ~ 10 6 . A Schwarzschild black hole 
is clearly a very poor oscillator. 



COMPLEX ANGULAR MOMENTUM 
APPROACH 



D. Useful approximations 

The quasinormal modes provide (at least in principle) 
a unique way of identifying black holes and deducing their 
mass and rate of rotation. Given this, it instructive to 
have simple approximations of the most important mode- 
frequencies. We can readily arrive at such expressions by 
recalling that the black hole problem is essentially one of 
scattering off a single potential peak (close to r = 3M) . 
It is commonly accepted that scattering resonances (the 
quantum analogues of quasinormal modes) arise for en- 
ergies close to the top of a potential barrier. This imme- 
diately leads to the approximation 



Re luq 



1 



3V3M 



1 + 



1 



(4.15) 



This is a good approximation of the fundamental 
Schwarzschild quasinormal (gravitational wave) mode for 
large I. For the imaginary part of the frequency — in 
quantum language: the lifetime of the resonance — the 
curvature of the potential at the peak contains the rel- 
evant information. Schutz and Will |llj used the WKB 
approximation to infer that 



Im ujq 



V3_ 
18M 



(4.16) 



As early as 1972 there was an intriguing suggestion 
(due to Goebel(l^]), that the black-hole resonant modes 
could have the following physical interpretation: a stand- 
ing wave could establish itself along the stable photon 
orbit at or near radius r — 3M. Of course this standing 
wave is not stable but would decay by radiating away 
energy. Unlike the quasi-normal modes, Goebel's stand- 
ing waves correspond to poles of the Green function in 
the complex angular momentum plane. Thus we need to 
extend our previous analysis to allow for complex values 
of I. That this leads to a powerful description of many 
scattering problems is well known |Lj ] . The complex an- 
gular momentum paradigm has been much investigated 
in acoustical and electromagnetic scattering, but has re- 
ceived scant attention in the context of black holes. 



A. Cross sections 

Let us begin by reviewing the theory of complex 
angular-momentum (CAM) scattering. Let F(l + 1/2) 
be any function which is analytic in the neighbourhood 
of the positive real / axis. Then, by Cauchy's theorem, 
we may write: 



E<-'>*w-i/ 



cos7rA 



dX 
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where A = I + 1/2. We apply this transformation to 
equation ( |3.16| ), writing 



2u 



c 



cos7rA 



dX. 



We deform the contour C away from the real I axis, 
rewriting / as the sum over the poles in Si and a back- 
ground integral (see Figure 15). 



f(0) = fp + 1b = 



-ITT 



E 



+ 



X n r n 



uj ' — ' cos A„ 

n 

l_ f X[Si-l]Pi(-cost 
2lo 



cos 7rA 



dX. (5.1) 



Here r n is the residue associated with the pole l n , defined 
by 

(5.2) 



I In 



in the vicinity of the nth pole. Even though we will not 
consider the 'background integral' further in this article, 
it is worth mentioning that it can be approximated using 
the saddle point method. 



FIG. 15: Integration contours in the complex A(= I + 1/2)- 
plane used in the derivation of the CAM representation for 
the scattering amplitude, equation (pM|). C is the original 
contour used in the integral representation for the scattering 
amplitude. The relevant contour for the background integral 
in the CAM picture is F. The Regge poles A n are all situated 
in the first quadrant. Their contribution is accounted for by 
the residue-theorem. 
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The poles of 5(A) are known as Regge poles and for 
the black hole potential one can show (Andersson and 
Thylwe, 1994) that they are all located in the first quad- 
rant of the complex A-plane. If If we compare the solution 
for one of the Regge poles to that of a quasinormal mode 
we see that they are rather similar. Both solutions sat- 
isfy purely "outgoing" wave boundary conditions. Hence, 
methods used for finding the quasi-normal mode frequen- 
cies can be adapted to find the Regge poles. One can 
show that these poles are all located in the first quad- 
rant of the complex A-plane. Typical results for the first 
few poles are shown in Figure [l6|. 



FIG. 16: The trajectory of the Regge poles Z n (w) (n = 
1,2,3,4) followed in the complex I plane from luM = (on 
the line I = -1/2) to coM = 5. 
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VI. PHYSICAL INTERPRETATIONS 

An important aspect of the CAM description of scat- 
tering is that each Regge pole has a clear interpretation. 
To realize this two approximations for the Legendre func- 
tions are useful: 

For 9 not close to and |A„| >> 1 we can use the 
formula 



Pl n (- cos t 



1/2 



Jo(A„(tt-0)) 



(6.1) 



when evaluating the Regge-pole sum (|d]). This is espe- 
cially interesting since we know that black-hole cross sec- 
tions show a prominent glory in the backward direction, 
see Figure [n} It is commonly understood that glories 
are characterized by Bessel-function type oscillations, cf. 
(| 3.25 ). When numerical results for the Regge pole with 



the smallest imaginary part (for a given ui) is used in (6.1) 
we get a good approximation to the glory oscillations in 
the black-hole cross section. 

Alternatively, we can use an asymptotic approximation 
for the Bessel function in (6.1). We then get 



iA„(7T-9)-i7r/4 , -«A„(7r-9)+i7r/4 

PiA- cos 9) « I , (6.2) 

V 27T A„ sin 9 

for |A„ sinf?| — * oo. From this formula it follows that we 
may interpret each Regge state as a combination of two 
surface waves travelling around the scattering centre (the 
black hole) in opposite directions. The angular velocity 
of each wave is proportional to 1/Re A n , and as they 
propagate around the black hole the waves decay expo- 
nentially. The imaginary part of A„ is clearly associated 
with the inverse of the "angular life" of each surface wave. 
It also follows, since Im A n > 0, that the amplitude of the 



first term in (6.2) is, in general, smaller that the second 
term. Only for 9 ~ ir do the two amplitudes have similar 
magnitude. Hence, one would expect interference effects 
to be more pronounced in the backward direction. More- 
over, it is easy to show that the anticipated diffraction 
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oscillations will have a period of 7r/Re A„. Again, this 
result approximates the features seen in the black-hole 
cross section rather well. 



TABLE I: "Angular life" and impact radius (R n ) for the first 
few Regge poles for uoM — 2. It should be remembered that 
the impact parameter associated with the unstable photon 
orbit at r = 3M is roughly 5.196M. 

n Impact radius (M) Angular life (degrees) 



5.194 
5.207 
5.234 



114 

38 
23 



We can also use the standard localization principle (cf. 
the definition of the impact parameter) 



Re A n ~ LoR n 



(6.3) 



to infer that the real part of each pole position is asso- 
ciated with the distance from the black hole at which 
the angular decay occurs. In the case of a Schwarzschild 
black hole one would expect such surface waves to be lo- 
calized close to the unstable photon orbit at r = 3M [or, 
strictly speaking, the maximum of the effective poten- 
tial in (3.4)]. This would correspond to R n = 3y/3M w 
5.196M. As can be seen from Table | the first Regge pole 
for various frequencies leads to a value of R n that is close 
to the impact parameter for the photon orbit. 



FIG. 17: The differential cross section for uiM = 2.0 as 
obtained from the phase-shifts (solid line) is compared to the 
contribution from each of the first three Regge poles (dashed 
lines). 
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the black hole. At the same time the waves decay at 
a rate related to the imaginary part of the Regge pole 
position. We have also seen that the glory oscillations 
that arise for large scattering angles in the black-hole 
case are naturally described in the CAM representation. 
In the specific example presented (for uiM = 2.0) a single 
Regge pole accounts for all large-angle structure in the 
scattering cross section. 



VII. QUANTUM EFFECTS 



A. Sample results 

Although the physical interpretations discussed above 
are suggestive and agree well with the partial-wave re- 
sults for the cross section we should also compute the ac- 
tual cross sections before assessing the usefulness of the 
CAM approach to black-hole scattering. In Figure |l7] we 
show the contribution to the the pole-sum fp from the 
first three Regge poles for u>M = 2.0. The results are 
compared to the partial-wave cross section as computed 
from the partial-wave sum (Figure |ll|). To obtain the 
Regge-pole contributions we have used the asymptotic 
formula ( |6.2j ) for the Legendre functions. 

From the data shown in Figure |l7] we make two ob- 
servations: i) For large scattering angles (9 > 40°) the 
pole sum in (5.1) is dominated by the Regge pole with 
the smallest imaginary part. Each consecutive pole gives 
a contribution that is roughly two orders of magnitude 
smaller than that of the preceding pole. This means that 
only one Regge pole need be included in a reasonably ac- 
curate description of the black-hole glory, ii) For smaller 
angles (9 < 40°) the first three terms in the pole sum are 
of the same order of magnitude. This is consistent with 
the interpretation of 1 /Im A„ as the "angular life" , see 
Table |. 

In conclusion, we have shown that the Regge states 
can be interpreted as surface waves that travel around 



We will now turn to another extreme scale in physics, 
the level where the classical theory of general relativity 
breaks down and it must somehow be married to quan- 
tum theory. We will not attempt to describe the ongoing 
attempts to formulate theory of quantum gravity in de- 
tail. Rather, we want to point out that one can gain some 
insights into quantum gravity (at the semiclassical level) 
using concepts and techniques very similar to those we 
have discussed above. In fact, several interesting effects 
bear a great resemblance to various scattering scenarios. 
Since we will only skim the surface of a vast area of re- 
search, we refer the reader to the monograph by Frolov 
and Novikov [n4| for references and more details. 



A. Hawking radiation 

Consider Figure [l|| A cloud of pressureless dust un- 
dergoes collapse to form a black hole. The wiggling lines 
represent high-frequency waves associated with a mass- 
less field (scalar, neutrino, photon, graviton). We sup- 
pose that the space is initially indistinguishable from 
Minkowski space and that the field is in the vacuum 
state. Now consider the evolution of a high-frequency 
wave packet that is a part of the spectrum of (Minkowski) 
vacuum fluctuations. We follow the evolution of the 
packet which, after the collapse, is localised just outside 
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FIG. 18: The Hawking effect: when dust collapses to form a 
black hole, high-frequency vacuum fluctuations near the hori- 
zon lead to radiation which can escape to infinity. 




the black hole's horizon at r = 2M. 

The packet can be spectrally decomposed into a com- 
ponent which is ingoing and another which (eventually) 
can escape to infinity. Hawking |ll| showed that, when 
viewed from infinity, the spectrum of the escaping ra- 
diation is thermal with temperature Tjj = h/ '(ks&irM). 
That is, if we Fourier-decompose the radiation escaping 
to infinity and use the standard definition of the vacuum 
state, we find that the total number of particles in the 
mode (n) is 

^ ( „) = |T (n) | 2 /(e^ BT » ) Tl) 

Here (n) is a general index which represents all quantum 
numbers and the sign is taken for negative for integer and 
positive for half-integer spin fields. 

We see the connection with scattering theory when we 
write the expression for the total (scalar) luminosity of 
the black hole due to Hawking radiation: 

oo 

H f uj du> v-^ / , m m / n , 9 

£ = ^ J / eWff _ 1 E(^ + 1 )i^Mi 2 - 

l =° 

There are similar formulas for quantum fields of other 
spins. For a stellar-mass black hole, this is 

Scalar = 7-44 x lO^h/M 2 

Reverting momentarily to cgs units, we can express the 
luminosity for a field of spin s as 

L s = a s ■ 4.09 x 1CT 17 ( ergs/sec. 

V M BH/ 

where Mq is the mass of the sun, and the a s are given 
in Table y. Clearly the Hawking effect is too small to 
ever detect for a black hole of near stellar mass. But the 
effect may be appreciable for much smaller black holes 



conceivably created in the early universe. Such primor- 
dial black holes would radiate at an appreciable level, and 
there have been speculations that one might be able to 
observe bursts of radiation as these black holes evaporate. 
As yet this effect has, however, not been observed. 



field 


spin s 


a a = LM 2 /h 


scalar 





7.44 x 10" 5 


neutrino 


1/2 


8.18 x 10" 5 


photon 


1 


3.37 x 1(T 5 


graviton 


2 


0.38 x 1(T 5 



TABLE II: Black hole luminosity as a function of field spin s 



B. Stress-energy tensors 

It must be emphasised that the Hawking radiation 
emitted by the black hole is not produced by the process 
of collapse but is feature of the field vacuum state in the 
new spacetime, which now contains a black hole. One of 
the lessons of quantum field theory in curved spacetime 
is that we must mistrust the abstract concept of 'parti- 
cle' and concentrate on physically measurable quantities 
such as the energy content of the field in a particular 
state, that is, the stress-energy tensor of the quantum 
field, written (T^"). It turns out that one can determine 
the stress-energy tensor via a differential operator acting 
on the black-hole perturbation Green's function that we 
have introduced earlier. 

It will be useful for what follows to define the Tolman 
local temperature 

Tloc - T H g 00 1/2 = T H I y/l - 2M/r. 
A massless scalar field at temperature T has 

(%"\ = ^T 4 diag(-3, 1,1,1) 

in flat space. A gas of scalar particles in equilibrium 
exterior to a spherical body should have 

(V) = ^locdiag(-3, 1,1,1) 

that is, T is replaced by the local temperature. This 
expression diverges as r — > 2M. 

The most obvious choice for the vacuum state exterior 
to a black hole is the state that, at infinity, resembles 
as much as possible the Minkowski vacuum state. This 
state, called the Boulware state, is the vacuum state ap- 
propriate for a field outside of a spherical body which is 
larger than its Schwarzschild radius, such as a neutron 
star. When (T^") is calculated in this state it has the 
asymptotic form 

~ -S T loc dia g(- 3 > !. M) r ^ 2M - 
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In other words, the energy density is negative infinite on 
H + and H , cf. Figure ^[ Note that this infinity has 
nothing to do with renormalisation: the divergences in 
(T^ 11 ) have already been removed using covariant, state- 
independent methods. The presence of a singularity in 
the stress-energy tensor, even when measured by a freely- 
falling observer, is unphysical and we conclude that the 
Boulware state cannot be the ground state of a field ex- 
terior to a black hole. 

If we are interested in a black hole formed by stellar 
collapse, we would like the stress-energy tensor to (at 
least) be regular on H + and to tend to the Minkowski 
vacuum tensor on I . The price we pay for this require- 
ment is that the tensor is singular on H . (Of course, 
in the case of stellar collapse, H~ does not exist.) In 
addition, at I + we have an outflow of thermal radiation 
as described in the preceding section. This state is called 
the Unruh state. 

If we demand that the stress-energy tensor be regular 
on both H and H + the Hartlc-Hawking state, and is 
appropriate for a spacetime with an eternal black hole 
in equilibrium with a thermal bath of radiation at the 
Hawking temperature. The stress tensor has the form: 



<v>* 



-2ft c diag(-3, 1,1,1) 
+finitc 



/2AA 6 / 6M 



\ r J 



black holes to spin very fast. Then the resultant space- 
time metric will no longer be spherically symmetric. The 
corresponding solution to Einstein's equations was dis- 
covered by Kerr, and the metric is usually written 



ds 2 



where 



r [dt — a sin 2 



sin 2 6 



[(r 2 +d 2 )^-adf 



4dr 2 
A 



2Mr + a 2 



p 2 d9 2 (8.1) 



p 2 = r 2 + a 2 cos 2 i 



Here a is a parameter representing the angular momen- 
tum per unit mass. When a equals zero the Kerr metric 
reduces to the Schwarzschild metric. The event horizon 
of a rotating black hole corresponds to the outer solution 
to A = 0, and is given by 



r+ = M + \/M 2 - a? < 2M. 



(8.2) 



In studying scattering from rotating black holes all the 
concepts we have introduced for non-rotating black holes 
remain useful. Thus we only need to comment on how 
these results are affected by the black holes rotation. It 
ig natural to begin by discussing the nature of light tra- 
jectories in the Kerr spacetime. 



where 'finite' is a correction everywhere finite and of or- 
der r~ 6 for large r. Despite its appearence, (T^) H is 
finite at r = 2M. 

We see then that the three 'vacuum' states correspond 
to different physical situations, and that it is not possi- 
ble to define a state which has the properties which we 
naturally associate with the vacuum. Although we re- 
jected the Boulware state because of infinite energy on 
the horizon, the Hartle-Hawking state also has infinite 
energy, which is contained in the heat bath at infinity. 
The Unruh state, the most 'realistic' state, is appropri- 
ate for an eternally evaporating black hole of constant 
mass. A real black hole will lose mass by evaporation, 
evolving to an unknown final state. This unsatisfactory 
situation will only be resolved by a full quantum theory 
of gravity. 



VIII. THE KERR BLACK HOLE 

So far we have only considered the simplest class of 
black holes, namely those described by the spherically 
symmetric Schwarzschild solution to Einstein's equa- 
tions. There are other kinds of black holes as well. In 
order to generalise our discussion to the case of greatest 
physical interest we must allow our black hole to rotate. 
Due to conservation of angular momentum during the 
gravitational collapse one would expect most newly born 



A. Null geodesies in the Kerr geometry 

A description of the general trajectories of a parti- 
cle moving in the Kerr geometry is considerably more 
complicated than the Schwarzschild case, but because of 
the axial symmetry one would still expect pe — for 
motion in the equatorial plane (p^ represents the four- 
momentum of a photon). If a particle is initially mov- 
ing in the equatorial plane, it should remain there. In 
the Schwarzschild case we could always, because of the 
spherical symmetry, orient our coordinate system in such 
a way that a study of equatorial trajectories covered all 
possible cases. For rotating black holes, we no longer do 
this and an equatorial trajectory is a very special case. 
Nevertheless, they provide a useful starting point for an 
exploration of particle motion around a rotating black 
hole. 

Because the Kerr geometry is stationary and axisym- 
metric, we will still have the two constants of motion 
p t = —E and p v = L z , the "energy measured at infinity" 
and the component of the angular momentum parallcll 
to the symmetry axis of the spacetime. Given this we 
can immediately deduce two equations of motion 

, dt (r 2 + a 2 ) 2 - a 2 A sin 2 6 2aMr , s 
P = ~d\ = SA E ~^A L ^ 



dip 2aMr A - a 2 sin 2 9 
P v = — = — E H — — : „ „ L 



d\ EA 



£Asin 2 6> z ' 



(8.4) 
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where S 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 6. 

We will restrict our attention to photons moving in the 
equatorial plane. Then the above formulae together with 
Pfj,p^ = and pg — lead to an equation for the radial 
motion that can be factorized as 



(r 2 



,2\2 



where we have defined 



a 2 A 



{E-V+)(E-V-), (8.5) 



2aMr±r 2 A 1 / 2 
V±{r) = — — —L, 



(r 



j 2 ) 2 -a 2 A J 



(8.6) 



By expanding these potentials in inverse powers of r 
we see that they fall off as 1/r as r — > oo (as in the 
Schwarzschild case), but the effect of rotation enters at 
order r~ 3 . This means that the rotation of the black hole 
has little effect on a distant photon. But as the photon 
approaches the black hole the potentials have a much 
stronger influence and we can distinguish two different 
cases. The way that the rotation of the black hole affects 
an incoming photon depends on direction of L z relative 
to the sense of rotation of the black hole. 

In the case when aL z > 0, when the photon moves 
around the black hole in a prograde orbit, we get the 
s itua tion illustrated in Figure 19. Then we can see from 
flgiel) that 



V- = at r = 2M 
V + = V-= " L 



•7) 



2Mr 4 



lo+Lz at r 



where we have defined the angular velocity of the horizon 

LU+. 



FIG. 19: The effective potentials for a photon moving in the 
equatorial plane of a rotating black hole. The figure illustrates 
the case when the photon has angular momentum directed in 
the same sense as the rotation of the hole (a = 0.5M). The 
corresponding figure for a retrograde photon (with aL z < 0) 
is obtained by turning this figure upside-down. 
03 
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Since the lefthand side of (3.5) must be positive (or 
zero) we can infer that the photon must either move in 
the region E > V + or in E < V_. In the first case the 
result is familar. An incoming photon from infinity can 



either be scattered by, or plunge into, the black hole. 
But what about the region V_ > E > that would also 
seem to be accessible (cf. Figure 19)? An analysis of this 



possibility requires some care. It turns out that it is not 
sufficient to require that E > 0, as one might intuitively 
think. The reason for this is quite easy to understand: E 
is the energy measured at infinity, and as we get closer to 
the black hole it becomes a less useful measure of what is 
going on. In order to understand the properties of light 
trajectories in the vicinity of a spinning black hole, we 
need an observer located close to the horizon to do our 
measurements for us. 

A convenient choice of local observer is one that has 
zero-angular momentum and resides at a fixed distance 
from the black hole (at constant r). Such Zero- Angular 
Momentum Observers (ZAMOs) were first introduced by 
Bardeen Q. (It should be noted that a ZAMO does 
not follow a geodesic, and consequently must maintain 
its position, say, by means of a rocket.) The suggested 
character of a ZAMO means that it must have a four- 
velocity 



= A 



: OjA, 



u r = u° = 



(8.9) 



The unknown coefficient A is specified by the requirement 



and we find that 



A 2 



-1 



9<p<p 



(8.10) 



(9ipt) 2 — gug vv 



.11) 



We are now equipped to address the question of pho- 
tons in the region V_ > E > in Figure |j| A ZAMO 
will measure the energy of a photon as 



-(ptu t +p v u*) = A{E-u;L z ) 



.12) 

This "locally measured" energy is the one that we must 
require to be positive definite, which means that we must 
have E > V + in Figure In other words, the possibility 
V- > E > is not physically acceptable and we can con- 
clude that the case aL z > only contains the same types 
of photon trajectories as we found in the Schwarzschild 
case. 

This is not, however, true for the case when aL z < 0, 
when the photon is inserted in a retrograde orbit around 
the black hole. (The potentials for this case are easily 
obtained by turning the ones in Figure [l9| upside down.) 
We then find from Q) that 



V+ = at r = 2M 



(8.13) 



and it is clear that some forward-going photons (that 
must lie above V+ according to our previous analysis) 
can have E < 0. That is, negative energy (as measured 
at infinity) photons may exist close to the black hole (for 
r < 2M in the equatorial plane)! As can be inferred 
from Figure [19 these negative energy photons can never 
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escape to infinity, but the fact that they can exist has an 
interesting consequence. 

Let us suppose that a pair of photons, the total energy 
of which is zero, are created in the region r + < r < 2M. 
The positive energy photon can then escape to infinity, 
while the negative energy one must eventually be swal- 
lowed by the hole. The net effect of this would be that 
rotational energy is carried away from the black hole, and 
it will slow down. This energy extraction process, that 
was first suggested by Penrose 0, can be extended to 
other objects. One can simply assume that a body breaks 
up into two or more pieces. If one of them is injected into 
a negative energy orbit the sum of the total energy of the 
remaining pieces must be greater than the total energy 
of the original body, since E is a conserved quantity. As 
in the case of photons, the extra energy is mined from 
the rotation of the black hole. But however exciting the 
possibility may seem, the Penrose process is unlikely to 
play a role in an astrophysical setting. 



B. The ergosphere 

As we have seen, there is a region close to a rotating 
black hole (r < 2M in the equatorial plane) where energy 
becomes 'peculiar'. This is the so-called ergosphere, and 
since there are many interesting effects (like the Penrose 
process) associated with it, it is worthwhile to discuss it 
in more detail. 

Let us consider a photon emitted at some r in the 
equatorial plane (9 = 7r/2) of a Kerr black hole. Assume 
that the photon is initially moving in the ±ip direction. 
That is, it is inserted in an orbit that is tangent to a 
circle of constant r. In this situation it is clear that only 
dt and d(p will be nonzero, and we find from ds 2 = that 




(8.14) 



From this we can see that something interesting happens 
if gu changes sign. At a point where g u = we have the 
two solutions 



dt Qtp^p 



and 



dt 



(8.15) 



In the Kerr geometry the first case corresponds to a pho- 
ton moving in the direction of the rotation of the black 
hole. The second solution, however, indicates that a pho- 
ton sent "backwards" does not (initially) move at all! 
The dragging of inertial frames has become so strong 
that the photon cannot move in the direction opposite 
to the rotation. Consequently, all particles must rotate 
with the hole, and no observers can remain at rest (at 
constant r, 9, (p) in the ergosphere. 

As the above example indicates, the boundary of the 
ergosphere follows from g u — 0. In the Kerr case we find 
that this corresponds to 



go = M ± \JM 2 - a 2 cos 2 6 



.17) 



From this follows that the ergosphere always lies outside 
the event horizon (even though it touches the horizon at 
the poles). 



C. Teukolsky's equation 

We want to extend our study of various scattering sce- 
narios to the Kerr case. To do so, we need to discuss 
perturbations of the Kerr geometry. It turns out that 
this problem is considerably more complicated than the 
Schwarzschild one. For example, the direct derivation 
of the equations governing perturbations of Kerr space- 
times by considering perturbations of the metric fails. 
It leads to gauge-dependent, and rather messy, formula- 
tions i n w hich one cannot readily separate the variables 
as in ( |3.3|) . 

A theoretically attractive alternative is to examine cur- 
vature perturbations. Using the Newman-Penrose for- 
malism, Teukolsky (1973) derived a master equation gov- 
erning not only gravitational perturbations (spin weight 
s = ±2) but scalar (s = 0), two-component neutrino 
(s = ±1/2) and electromagnetic (s = ±1) fields as well. 
In Boycr-Lindquist coordinates and with the use of the 
Kinnersley null tetrad, this master evolution equation 
reads 



(r 2 +a 2 ) 



2\2 



a 2 sin 2 9 



_ _ 4Mar_ 
dtt'f t— d t ^ 



-2s 



M(r 2 - a 2 ) . 
r : h la cos ( 



d f * 



A~ s d r (A e+1 d r f) + —-de (sinflfytf) 
v ' sin 9 



sin 

2 + 2 



v s+l 

a* 
A 



a(r — M) i cos ( 



A 



s 2 cot 2 9-s)* = 



(8.18) 



A - a 2 sin 2 9 = 



(8.16) 



The actual meaning of W for various spin-fields is rather 
complex, so we will only worry about two special cases 
here. Firstly, for s — \& represents the scalar field (<&) 
itself and in the limit a — > we recover the Schwarzschild 
scalar wave equation. For s — ±2 the * corresponds 
to the Weyl curvature scalars and ^4 that directly 
represent the gravitational-wave degrees of freedom. 

The great brea kthrough that followed Teukolsky's 
derivation of (3.18) was that one could now separate the 
variables also for Kerr perturbations. For our present 
purposes it is sufficient to note that this essentially cor- 
responds to assuming that i) the time-dependence of the 
perturbation is accounted for via Fourier transformation, 
and ii) there exists a suitable set of angular function that 
can be used to separate the coordinates r and 9. In the 
case of scalar perturbations, the angular functions turn 
out to be standard spheroidal wavefunctions. Knowing 
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this we assume a representation (for each given integer 

771 ) 



$ = / doj 



1=0 



Ri m (r,uj)Si m (9,uj) 



(8.19) 



where it should be noted that the angular functions de- 
pend explicitly on the frequency to. That is, they are 
intrinsically time-dependent functions. After separation 
of variables, the problem is reduced to a single ordinary 
differential equation for Ri m (r,uj). This equation can be 
written as 



d Rl n 

dri 



K 2 + (2amuj - a 2 u? - E)A dG q2 
(r 2 +a 2 ) 2 



dr 

(8.20) 

where K = (r 2 + a 2 )uj — am, G — rA/(r 2 + a 2 ) 2 , 
and the tortoise coordinate is defined from dr* — 
(r 2 + a 2 )/ A dr. The variable E is the angular separa- 
tion constant. In the limiting case a — > 0, it reduces to 
1(1+1), and for nonzero a it can be obtained from a power 
series in auj. It should be noted that E is real valued for 
real frequencies. 



D. Quasinormal modes 



The physical solution to (8.2C ) is defined by the asymp- 
totic behaviour (cf. the Schwarzschild result) 



Rlr. 



as r 
as r 



+oo . 



(8.21) 

where u+ = a/2Mr+ is the angular velocity of the event 
horizon. 

From this we see that we can define the quasinormal 
modes of a Kerr black hole in the same way as we did 
in the spherically symmetric case. Furthermore, these 
modes can also be calculated using Leaver's continued 
fraction method. The results can be summarised as fol- 
lows: Recall that in the non-rotating case the modes 
occur in complex-frequency pairs u q and — uj q (the bar 
denotes complex conjugation). This is apparent in Fig- 
ure |l4|. As the black hole spins up, each Schwarzschild 
mode splits into a multiplet of 21 + 1 distinct modes (in 
analogy with the Zeeman splitting in quantum mechan- 
ics) . These modes are associated with the various values 
of to, where — I < m < I, which determine the modal 
dependence on the azimuthal angle through e imv . As is 
straightforward to deduce, modes for which Re u q and m 
have the same sign are co-rotating with the black hole. 
Similarly, modes such that Re u q and m have opposite 
signs are counter-rotating. The effect that rotation has 
on the mode-frequencies can, to some extent, be deduced 
from this fact. 

Let us first consider the counter-rotating modes: These 
will appear to be slowed down by inertial frame dragging 



close to the black hole. Hence, their oscillation frequen- 
cies will tend to decrease as a — > M. At the same time, 
numerical calculations show that the damping rate stays 
almost constant. For the co-rotating modes, the effect is 
the opposite. Frame-dragging tends to increase the fre- 
quencies. Additionally, the modes become much longer 
lived. The available numerical results for co-rotating 
modes are well approximated by (cf. fTpf ) 



Re wo 



1 

M 



1 



and 



Rim — 



Im ujq 



(1 - a/M) 
AM 



63 

Too 



9/10 



(1-a/M) 3 / 10 



.22) 



100 ^ 



a/M) 



3/10 



.23) 

From this we can see that the mode becomes undamped 
in the limit a — M. One can actually show (since the 
case a = M is amenable to analytic methods) that there 
exists an infinite sequence of real resonant frequencies 
with a common limiting point in that case. The limiting 
frequency is u> — m/2M, which we will later show to be 
the upper limit for so-called super-radiance. That the 
modes become undamped can be understood from the 
fact that the angular frequency of an extreme Kerr black 
hole is 1/2M. As a — > M the long lived quasinormal 
modes essentially rotate uniformly with the black hole, 
and as a consequence they do not radiate strongly. 

Provided that the modes of close to extreme rotating 
black holes will be excited by some realistic astrophysical 
process, the fact that these modes can be very long-lived 
would greatly improve the chances for detection with fu- 
ture gravitational- wave detectors. Hence, it is of interest 
to investigate the excitation of these long lived modes. 
Recent work has shown that the modes tend to be harder 
to excite than the short-lived Schwarzschild modes, but 
that the slowly damped modes nevertheless dominate the 
emerging signals. An example of this is shown in Fig- 
ure M 



E. Superradiant scattering 



Given the prescribed asymptotic behaviour (8.21), to- 
gether with that for the complex conjugate of Ri m and 
the fact that two linearly independent solutions to ( 8.20| ) 
must lead to a constant Wronskian, it is not difficult to 
show that 



(l-muj + /uj)\T\ 2 = l-|Sf 



(8.24) 



where we have introduced the transmission and reflection 
coefficients as in the Schwarschild case; 



ITI 



1 



From this result, it is evident that the scattered waves 
are amplified (l^ 2 > 1) if 

lu < muj + (8.25) 
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FIG. 20: The response of a near extreme Kerr black hole af- 
ter a Gaussian scalar wave pulse has impinged upon it. The 
main features are the same as in the Schwarzschild case (Fig- 
ure |l3| ) , but here the quasinormal mode ringing is much slower 
damped. 
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This is known as superradiance and it is the wave- 
analogue to the Penrose process that we described earlier. 
Its existence implies that it would in principle be possible 
to mine a rotating black hole for some of its rotational 
energy. 

In Figure |2l], we show a sample of results for the re- 
flection coefficient in the case when I = m = 2. These 
r esult s were obtained by a straightforward integration of 
( |S.20 ) and subsequent extraction of 1Z. The maximum 
amplification in this case is a minute 0.2 %. Similar re- 
sults for electromagnetic waves and gravitational pertur- 
bations show that the maximum amplification is 0.3% for 
scalar waves, 4.4% for electromagnetic waves and as large 
as 138% for gravitational waves. 

Given the results in Figure ^], it is worth pointing out 
that they agree with standard conclusions regarding the 
apparent "size" of a rotating black hole as seen by differ- 
ent observers. A rotating black hole will appear larger to 
a particle moving around it in a retrograde orbit than to a 
particle in a prograde orbit. This is illustrated by the fact 
that the unstable circular photon orbit (at r — 3M in the 
non-rotating case) is located at r = iM for a retrograde 
photon, while it lies at r = M for a prograde photon. 
The results in Figure 2l] illustrate the same effect: In our 
case, we have prograde motion when to/m is positive and 
retrograde motion when uj/m is negative. The data in 
Figure ||l| correspond to m — 2, and the enhanced reflec- 
tion for positive frequencies as a — > M has the effect that 
the black hole "looks smaller" to such waves. Conversely, 
the slightly decreased reflection for negative frequencies 
leads to the black hole appearing "larger" as a — > M. A 



sample of results showing this effect is in Table 111 . 



F. Scattering of waves by Kerr black holes 

One can analyse the scattering of monochromatic 
waves by a Kerr black hole in much the same way as 



FIG. 21: Reflection coefficient for different values of a in the 
case I = m = 2. Upper panel: As a increases, there is a clearly 
enhanced reflection of prograde waves {ui > in the figure) 
while the overall reflection of retrograde waves {ui < in the 
figure) decreases somewhat. Lower panel: A blow-up of the 
result for prograde waves unveils a maximum amplification 
due to superradiance of 0.187%. 
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TABLE III: The apparent size 6 of a Kerr black hole as viewed 
along the rotation axis. The values are all for a = 0.9M 
and are estimated from the total absorbtion cross section 
<r abs ~ irb 2 . Positive frequencies co-rotate with the black 
hole whereas negative ones are counter-rotating. The values 
should be compared to b — 5.2M for a Schwarzschild black 
hole. 
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we approached the Schwarzschild case. The same is true 
also for wave fields other than scalar waves. The analysis 
of electromagnetic and gravitational wave scattering is, 
in principle, identical to that for scalar waves. However, 
in the case of gravitational waves an additional complica- 
tion enters: Gravitational waves come with two different 
polarisations. This means that the scattering amplitude 
consists of a sum of their individual contributions, and 
that the cross section may show features due to interfer- 
ence between these two contributing terms. This effect 
has not been explored in detail as yet. 

Similarly, despite a few studies, the full details of scat- 
tering from rotating black holes remain to be understood. 
Scattering of gravitational waves incident along the axis 
of symmetry of a Kerr black hole show that the scatter- 
ing cross section depends in a complicated way on the 
rotation parameter a. One would essentially expect Kerr 
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scattering to be different because of two effects that do 
not exist for nonrotating holes: superradiance and the 
polarization of the incident wave. For incidence along 
the symmetry axis of a rotating black hole one can have 
either co- or counter-rotating waves. The two cases lead 
to quite different results. Although the general features 
of the corresponding cross sections are similar they show 
different structure in the backward direction. This pos- 
sibly arises because of the phase-difference between the 
two polarisations of gravitational waves. As for superra- 
diance, it has been suggested that it tends to wash out 



interference minima. Further details can be found in the 
book by Futterman, Handler and Matzner (l9||. 
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